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A recently introduced capillary-wave description of binary-alloy solidification is generalized to 
include the procedure of directional solidification. For a class of model systems a universal dispersion 
relation of the unstable eigenmodes of a planar steady-state solidification front is derived, which 
readjusts previously known stability considerations. We, moreover, establish a differential equation 
for oscillatory motions of a planar interface that offers a limit-cycle scenario for the formation of 
solute bands, and, taking into account the MuUins-Sekerka instability, of banded structures. 



PACS numbers: 68.35.Dv,81.10.Aj,05.70.Np 

INTRODUCTION 

The main feature of a recently introduced capillary- 
wave model [I] for the solidification of a dilute binary 
alloy is the use of the interface position as a basic field 
variable, in addition to the concentration of the solute 
component. In the present paper this approach will be 
generalized to cover also the description of directional 
solidification, especially with regard to the rapid-growth 
regime. As outlined in reviews by Langer [2] and by 
Miiller-Krumbhaar et al. |3], in directional solidification 
the growth of a crystal is accomplished by pulling it in 
opposite direction of an externally applied temperature 
gradient. We will mainly consider the case of a constant 
temperature gradient, which enters via a driving force in 
the equation of motion for the interface position. This 
form of description arises in the limit of an infinite heat 
conductivity from a more general model, involving energy 
density as an additional field variable. In general, such a 
model would allow to include the effect of heat diffusion. 

As a first application of our approach we scrutinize the 
possibility of stationary motions of a planar solidification 
front. The stability of such a front has been investigated 
in the rapid-growth regime with increasing regard of non- 
equilibrium effects by MuUins and Sekerka HI , Coriell and 
Sekerka ^, , and by Merchant and Davis ,6J. In Refs. [5 
and [6] new oscillatory interface instabilities have been 
discovered in addition to the previously-known MuUins- 
Sekerka instability. Our own approach demonstrates that 
these effects are closely related to an instability, found 
by Cahn ff^ in grain-boundary motion. The threshold of 
this instability represents a border line between regimes 
of steady-state and of non-steady-statc motions of the 
solidification front. 



Non-steady interface motions operate in generating the 
periodic growth of layers with alternating homogeneous 
and dendritic micro-structures in binary alloys. This so- 
called banded structure occurs in many metallic alloys, 
as described in the review [S] by Garrard et al. who also 
offer a phenomenological explanation of the effect, using 
a quasi-stationary approximation. In a more microscopic 
treatment, Karma and Sarkissian 9J pointed out that the 
banding phenomenon is due to relaxation oscillations of 
the solidification front. This behavior, described in more 
detail in Ref. [10] . was derived from numerical solutions 
of the diffusion equations for the solute concentration in 
a dilute binary alloy and for the temperature, supported 
by non-equilibrium boundary conditions, as formulated 
by Aziz and Bocttinger Starting from a phase-field 
model, Conti has performed one- and two-dimensional 
simulations, describing the generation of solute bands 
fl2, , and of banded structures , respectively. He also 
has confirmed in Ref. [T3| the observation by Karma and 
Sarkissian that the inclusion of heat diffusion leads to 
an increasing suppression of the formation of bands with 
decreasing heat conductivity. 

In the application of our capillary- wave description we 
are going to analyze the simplest-possible model, which 
shows the banding effect. We, accordingly, consider the 
dynamics of a planar interface, neglecting heat diffusion, 
and assuming an overall constant diffusion coefficient for 
the solute component. For a class of model systems with 
arbitrary equilibrium profiles of the solute concentration 
these properties lead to an integro-differential equation 
for the interface position. In case of a sufficiently small 
temperature gradient, realized in many experiments, this 
equation can be reduced to the differential equation of a 
damped nonlinear oscillator [15 . Stationary solutions of 
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this equation turn out to exist in some region, limited 
by the threshold of the Cahn instability. This instability 
is attended by an oscillation, blowing up in the unstable 
regime to a limit-cycle behavior of the solidification front 
and of the solute concentration at the interface. Close 
to the stability threshold the transition from uniform to 
periodic solutions can analytically be evaluated by the 
Bogoliubov-Mitropolsky method \16\. The nature of the 
periodic solutions can be tuned from almost- harmonic to 
distinctive relaxation oscillations by changing the pulling 
velocity from the stability threshold to values deep inside 
the unstable regime. 

A benefit of our approach is that, apart from solving 
the oscillator equation in the unstable regime, all steps of 
the procedure could be accomplished analytically, which 
deepens our understanding of the banding effect. The 
discussion of micro-segregation effects at an oscillating 
solidification front requires to consider the stability of 
such a front in the transverse direction. Since, however, 
our results reveal an almost stationary behavior of the 
interface motion in the so far barely understood low- 
velocity regime [lOj . we presently only complement the 
established limit-cycle scenario by the standard MuUins- 
Sekerka procedure. Then, in some window of the model 
parameters, the high- and low-velocity sections of a limit 
cycle are located inside the Mullins-Sekerka stable and 
unstable regimes. As a result, a dendritic microstructure 
will develop in the low- velocity bands, which we consider 
as a kind of noise on the more macroscopic scale of the 
periodic array of the widely flat bands. 

CAPILLARY- WAVE MODEL 

The effective Hamiltonian of our capillary-wave model 
is a functional of the interface position Z(x, t) and of 
the excess concentration C{r,t) of the solute relative to 
its value Cs in the solid phase. In terms of these field 
variables the effective Hamiltonian has the form 



H ^ - I (fx {dZf 



C-U{z-Z) 



(1) 



established already in Ref. \X^. It determines all static 
properties of the system in thermal equilibrium at some 
fixed temperature Ts < Tm where Tm denotes the 
melting temperature of the solvent, showing up in the 
temperature-concentration phase diagram. Fig. 1. The 
input quantities in the Hamiltonian ([I]) are the surface 
tension a, the coupling parameter 
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FIG. 1: Temperature-concentration phase diagram, showing 
the liquidus and solidus lines Tl{C),Ts{C) which meet at 
Tm- The values Cl and Cs refer to the temperature Ts- 



involving the solute concentration Cl in the liquid phase, 
the latent heat L per unit volume, the miscibility gap 



AC = Cl - Cs 



(3) 



visible in Fig. 1, and the solute-concentration profile in 
thermal equilibrium, 



U{z- Z) ^Ce{z-Z) 



(4) 



Whereas the expression ^ has been derived in Ref. [T], 
Eq. Q directly follows from the equilibrium condition 
5H/5C = 0. 

From Ref. [T] we also adopt the equations of motion 



dfZ = k \ F - 



SH 

Jz 



dtC 



1 5H 



K 5C 



(5) 



where the rate A measures the interface mobility, and 
D is the diffusion coefficient of the solute, here assumed 
to have an overall constant value. The externally applied 
temperature gradient S, and the pulling velocity Vp enter 
via the driving force 



F ^ L 



Ts-T 
Tm 



(6) 



where, in terms of the temperature Tp at the steady-state 
position Z{t) ~ Vp t, 



OCl 
dT 



L 1 
AC 



(2) 



T = Tp + S{Z - Zp) , Zp = Vpt 



(7) 



In Appendix A we will consider a more general model, 
which includes energy density as an additional field. We 
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also will show that the reduced model equations ([ij - Q 
emerge in the limit of an infinite heat conductivity. 

A dimensionless form of the model equations ([5| can 
be obtained by adopting from Ref. [1 the mappings 



1 D f 

-r^r , T^t^t , ^F^F 



(8) 



where, in the present context, the length ^ is defined by 
AC 1 



^ 2 [/'(O) 
In terms of the dimensionless quantities 



(9) 



= ^ f — \' 



Act 



eL s 



a Tm 

the resulting equations of motion read 



(10) 



-dtZ ^ F{Z)~-i / dzU'iz- Z)[C -U{z~ Z)] , 

dtC^DV^[C~U{z^Z)]) , (11) 
with the driving force F given by 



F^Fp-m\Z-Zp) , Fp = ^%^ . (12) 



<T T, 



As a first application of Eqs. (Ill and (12) we now will 



consider the steady-state motion of a planar interface 
with velocity Vp in z-direction. 



STATIONARY PLANAR GROWTH 



Measuring velocities in units of the diffusion velocity. 



(13) 



the stationary growth of a planar solidification front is 
described in the co-moving frame z = vpt + C by the 
equations 



1 



vp^Fp + Gpivp)~GpiO) , 



(14) 



/ + 00 
dCU'{C)Cp{C;vp) , 
-oo 



CpiCvp) 



dCU'iC) expbp(C'-C)] 



These equations are identical to those, derived in Ref. [T] 
for the case of solidification by under-cooling the liquid 



phase from Tg to Tp. In particular, the result (14 1 is 
independent of the parameter m^, which only enters in 
discussing the stability of the planar morphology. 
For perturbations of the form 

h{-K,t) = Z{x,t) ~vpt , (15) 
c(x, C, t) = C(x, C, t) - C7p(C; vp) + C'piQ- vp)hi^, t) , 

the resulting equations of motion read 



1 

-dth = {d^-m^)h- dCC/'(C)c(x,C,0 , (16) 

P J-oo 

dtc = VP d^c + {dl + d^)c+ [C'p{dt - d^) + U'd^] h 

where = V^ — These equations have eigensolutions 
of the form 



/i(x, t) = /i(q, w) exp (iq • X + a;t) , (17) 
c(C, X, t) = c(C, q, w) exp (iq • x + tot) , 

which, after elimination of the component c, lead to the 
eigenmode dispersion relation 



P 



+ + - vp[Gp{vp + A) - Gp{vp) 



[Gpivp + X) + GpiX)] (18) 



Vp + 2A 

where, deriving from the equation of motion for c. 



A = -{vp/2) + V(«p/2)2 +w + q2 



(19) 



The result (18 1 is similar to that, established in Ref. |T] 
and only differs by the additional term m^. 

Inspection of the low- q,u} behavior of the dispersion 



relation (18) leads to identify an eigenfrequency wi(g), 
which captures the MuUins-Sekerka instability. Since the 
parameter m acts as a long-wavelength cutoff, the wave- 
number threshold qc for this instability is shifted from 

= at m = to some finite value, determined by the 
relations u!i{qc) = a;5.(9c) = 0. Elimination of qc then 
generates the neutral stability curve of the instability in 
form of a function vp{"f), with a parametric dependence 
on m. An explicit form of this neutral line will later be 
derived for a specific expression of U{z — Z). 

A second branch a;2 (<?) comprises an instability, similar 
to that, discovered by Cahn in the process of grain- 
boundary motion. This instability is characterized by a 
gap at q = 0, determined by the relation 



^2(0) 

n{vp) 



2m^ 



± 



F'p{vp)n{vp) ^ [F'pivpMvpW 



4m4 



4 



n{vp) 



dvp 



Gpivp) + GpiO) 



2vp 



-1/2 



(20) 



In the limit to — >■ a single nonzero value of the gap 
survives, which is identical to that, found in Ref. [T]. 
For TO 7^ the neutral stability curve Fp{vp) = of the 
Cahn instability is attended by an oscillation of period 
U,(vp). Similar oscillatory instabilities have previously 
been observed by Coriell and Sekerka [5j, and later by 
Merchant and Davis [S]- However, the neutral line, found 
in Ref. [!)', differs from ours, which, we conjecture, arises 
from an unsettled generalization of the Gibbs-Thomson 
relation. We next will demonstrate that the instability. 



described by Eq. (20), acts as a seed for the limit-cycle 



behavior in the unstable regime. 



NON-STATIONARY PLANAR GROWTH 



For general unsteady motions of a planar interface the 



equations of motion (11 1 reduce to the form 



Z{t) ^Fp- m^[Z{t) - Zp{t)] 



+ 00 



dz U'(z - Z{t))[C{z, t) - U{z - Z{t))] 



{dt-dl)C{z,t) = -U"{z-Z{t)) . 



(21) 



We are mainly interested in the late-stage behavior of 
Z{t), and, therefore, are going to replace in the first of the 
equations (21 ) the solution C(z, t) of the second equation, 
complemented by the boundary condition C{z, — oo) = 0. 
This leads to the expression 



C{z,t) 



t r-\-oo 

dt' / dz' d^^g{z^ z\t-t') 



U'{z'-Z{t'}) , (22) 



involving the Green function 



r+°° 

Q{z,t)^ / — eyipi-k^t + ikz) . (23) 
J-oo 27r 

After the variable substitutions 

C = z- Z{t) , C' = ^'- Z{t') , (24) 
and expansion of Z(t') around Z{t), we obtain 



GiC -C + z{t)~zit'),t-t') 
dk 



(25) 



27r 



exp [-k\t - t') + ik{C - C') + - 



■ exp 



-rkY,^-^{t-trdr'v{t) 



using the notation 



v{t) = Z{t) ^vp + h{t) 



(26) 



For m? <C 1, the higher-order contributions in n are 
increasingly negligible, as seen from the scaling procedure 



h ^ m ^/i , dt ^ rn^dt 



(27) 



which leaves v{t) invariant, and attaches a factor to^"~^ 
to the contributions (x d"~^v{t). In the so-called quasi- 
stationary approximation all terms of order n > 2 are 
neglected. The scenario, developed by Garrard et al. [S], 
is based on this procedure, applied to a phenomenological 
model where a low-velocity dendritic branch is added to 
the curve F — Fp{v) for a planar interface. Without this 
additional dendritic branch all trajectories in the F,v - 
plane would inevitably run to t; = 0. The phase-field 
simulations by Conti [12] effectively include all n-order 



terms in Eq. (25 1, and for the planar interface lead to the 



appearance of limit cycles, which are well separated from 
the line v = 0. An almost identical behavior arises, if in 
Eq. (25 1 we just include the term n = 2 , thereby going 



one step beyond the quasi-stationary approximation. 

After expansion of the exponential in Eq. (25) up to 
n = 2, and collection of all terms in Eq. (22) depending 
on t', integration over t = t — t' yields 



ik j dr 
'o 

1 



1 — ik — i} 
2 



exp [(— fc^ + ik v) t] 



1 



k — iv (k + ie){k — iv)^ 



(28) 



The shift -I- ie in the denominator of the last term arises 
from including a term — ev in the preceding exponential, 
which regularizes the singular point A; = at the upper 
bound of the integral. 

Next, we take care of the, so far, ignored contribution 
i/c(C — C') in Eq. (25 ), and perform the integrations over 



k separately for the two final contributions in Eq. ( 28 1 
The resulting equations 



5 



+ 00 



^exp[^fc(C-C')]^ 

ZTT k ^ IV 



e(C-C')cxphz;(C-C')] , 



—00 
,2 



1 9 
1 92 1 



+00 



dk exp [zfc(C - CO] 
27r (/c + i£)(fc — iv) 



- {e(c - O exp hz;(c - C')] + 0(C' - 0} 

have, finally, to be multiplied with U'{C,') and integrated 



over C in order to evaluate the expression (22 1 for the 



solute concentration in the assumed approximation. 
In terms of the stationary concentration profile Cp, 



presented in Eqs. (14 1, the result for C{z,t) reads 



Ciz,t) = CpiC;v)+v^^^[CpiC;v)+Cp{C;0)] ■ (30) 



Insertion of this into the first equation in Eqs. (21 ) leads 



to the closed equation of motion for Z{t) in the form 



Fp - m?{Z - Zp) 



(31) 



+ Gp{v) - Gp(0) + v]-~[Gp{v) + Gp(0)] 

2 OV^ V 



where Gp{v) has been defined in Eqs. (14 1. For v = vp 
the result ([3l|) consistently reduces to the first line in 



Eqs. (14 1. Subtracting the latter from Eq. (31), we find 



for the displacement 



We mention that in the limit v{t) = vp + h(t) — ^ 
the inertial term shows, after application of the scaling 



(29) procedure ( 60 1 , the behavior 



M{h)h (X ^ 'h 



(35) 



The most singular terms in higher-order contributions 
turn out to carry a pre- factor {m^ /v^)"~^ , so that our 
oscillator equation is only valid for velocities above the 
cross-over line 



T 2 

V (X m 



(36) 



We, furthermore, observe that, for small /i, the second 



definition in Eq. ( 34 1 implies the behavior 



R{h) 



P 



G'p{vp) 



h + 0{h^ 



(37) 



for the friction term in Eq. (33), so that, due to the first 



line in Eqs. ( 14) 



M{0)h + F'p{vp)h + m^h + 0{h^) = , 
in agreement with our linear stability analysis. 



(38) 



The nonlinear differential equation ( 33 ) obviously has 



the trivial solution h{t) — 0, which, however, according 
to Eq. ( 20 ) , is unaffected by the Cahn instability in the 



regime Fp{vp) > 0. In the regime F'p{vp) < we will 
find solutions h{t), showing an oscillatory behavior in 
the limit t — 00. This behavior is shared by the solute 
concentration C{Z{t),t) at the oscillating interface, as 



can be seen from Eq. (30), taken at C = 0. 



h{t) EE Z{t) - Zp(t) 

the simpler differential equation 



(32) 



M{h{t)) h{t) + R{h{t)) + h{t) = (33) 
where we have introduced the mass and friction functions 



M{h) 



1 &^ 

2 dv^ 



Gp[vp + h) + Gp{Q) 
Vp + h 



1 ■ 



R{h) = -h-Gpivp + h) + Gp(vp) . (34) 
P 



Together with these definitions, Eq. (33) represents one 



of the central results of the present paper. It describes a 
damped nonlinear oscillator, general properties of which 
have been discussed in Ref. [TBI. 



LIMIT-CYCLE SOLUTIONS 



The equations ( 30 ) and ( 33 ) are valid for a whole class 



of models with varying equilibrium-concentration profiles 
U{z — Z). In order to obtain explicit solutions h(t) and 
G{Z{t),t), we choose the model 

U{z-Z) = Q{Z - z)exp{z - Z) (39) 
-I- e(z-Z)[2-exp(Z-z)] , 

derived in Ref. [T] from a two-parabola phase-field model. 
As also explained in Ref. [T] , the choice ( 39 ) leads to the 
expressions 



CpiO,v) 



1 



Gp{v) 



{v + iy 



(40) 



6 





FIG. 2; The functions M(h),R{h) for 7 = 0.01, p = 100, 
and vp = 0.3, resulting from the model (391 of the solute 
concentration. 



FIG. 3: Solutions h{t) for 7 = 0.01, p = 100, m ■ 
pulling velocities vp = 0.522, and v = 0.52. 



0.003, and 



The latter result allows us to determine the quantities 
( 34 1, which e. g. for 7 = 0.01, p = 100, wp = 0.3 have the 
form, shown in Fig. 2. A conspicuous property of the 
function M{h) is its monotonous growth with decreasing 
velocity, which decisively affects the solutions h(t) in this 
regime, and, therefore, supports the procedure to include 
the inertial term in the oscillator equation (33). Another 
implication of Fig. 2 is that, for the present choice of 
the model parameters, the function R{h) is negative in 
some finite region where the solution h{t) = is unstable. 
The first result in Eqs. (40), finally, permits to calculate 



solution h(t) of Eq. (33) has been found 



the solute concentration (30) at the interface, once the 



Numerically obtained solutions h{t) for the parameter 
values 7 — 0.01, p — 100, m = 0.003, and for the pulling 
velocities vp — 0.522 and vp ~ 0.52 are shown in Fig. 3. 
The threshold condition Fp{vc) = generally defines a 
critical velocity, which in the present case has the value 
vc ~ 0.521. Above vc the oscillating trajectories h{t) 
converge to the value h(oo) — whereas below vc they 
approach a limit cycle. Fig. 4 shows the same behavior 
further away from vq, so that, comparing these figures, 
one observes a kind of critical slowing down. 



In the regime \vp — vc\/vc ^ 1 the envelopes in 
Fig. 3 can be calculated analytically by the Bogoliubov- 
Mitropolsky procedure the application of which to 
the present case is described in Appendix B. The result 



for the solution of Eq. ( 33 ) then is found to be 



h{t) = a{t)costp{t) 



(41) 



where ipi'^) is a rapidly oscillating phase, and a{t) is an 
amplitude, obeying the differential equation 



da 
'dt 



= -pi a - p^a 



(42) 



Here, pi = ri{vp — vc), and the coefficients ri,p3 are 
determined by the values of the model parameters 7,p, to. 



Eq. ( 42 ) has the solution 



a{t) = ao 



1 



P3 2 
— Oo 

Pi 



exp [2pi t] — 



P3 2 
— ao 
Pi 



-1/2 



, (43) 



which for pi > and pi < describes the envelopes in 
Fig. 3. The asymptotic value of the limit-cycle amplitude 
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FIG. 5: Orbit of the cycle for 7 
and vp = 0.5. 



0.01, p = 100, m = 0.003, 



in Ref. [12]. From the associated behavior of C{Z{t),t) 
we, moreover, see that the transitions between high- and 
low-concentration layers are joined by large-acceleration 
sections. As already pointed out by Carrad et al. [5], 
this explains the appearance of relatively sharp interfaces 
between these layers. Fig. 8, finally, presents our result 
for the orbit of the limit cycle belonging to Fig. 7. 



FIG. 4: Solutions h{t) for 7 = 0.01, p = 100, m 
pulling velocities vp = 0.53, and vp = 0.51. 



0.003, and 



shows the critical behavior a(oo) = \/ ^Pi/ Ps- In the 
marginal case pi = Eq. ( 43 1 implies the algebraic decay 



ait) = 



(44) 



The rapid oscillations in a fully developed limit cycle, 
gleaming through in Fig. 4, are most suitably analyzed 
by numerical computations. A first result is the orbit of 
a limit cycle in the h,h - plane, shown in Fig. 5 for the 
parameter values 7 = 0.01, p = 100, to — 0.003, and the 
pulling velocity vp — 0.5. The related oscillations of the 
trajectories h{t),h{t), and of C{Z{t),t) are displayed in 
Fig. 6. Since the term m?h{t) measures the temperature 
at the oscillating interface, this quantity is effectively also 
included in Fig. 6. 

By reducing the pulling velocity to the value vp = 0.3 
at constant parameters ^,p,m, one obtains the shape 
of the trajectories h{t),h{t),C{Z{t),t) deeper inside the 
limit-cycle regime. The results for h{t), h{t), displayed in 
Fig. 7, are remarkably close to the findings by Conti 



BANDED-STRUCTURE FORMATION 

The layer formation, induced by the above limit-cycle 
solutions is unaffected by the Mullins-Sekerka instability, 
deriving from Eq. (18). This follows from Fig. 9, which 



shows the neutral stability lines, enclosing the unstable 
regions of the Cahn and the Mullins-Sekerka instabilities, 
and the projection of the limit cycle in Fig. 8. We have 
to point out, however, that the form of the the Mullins- 
Sekerka neutral line is only an approximate one, since it 
is related to a steady-state reference motion with velocity 
Vp. The approximation seems, however, to be acceptable 
due to the almost stationary behavior of h(t) in Fig. 7 
at low velocities. We, accordingly, expect that Fig. 7 
induces the formation of precipitation-free solute bands. 

The approximation for the Mullins-Sekerka neutral line 
is apparently more justified for the limit cycle, belonging 
to Fig. 10. In this case, the projection of the cycle enters 
the unstable region of the Mulins-Sekerka instability, as 
seen in Fig. 11. Accordingly, the interface will develop a 
dendritic microstructure at low velocities, which dissolves 
again in the high- velocity regime. This is just what one 
expects to happen in the creation of banded structures, 
and it is in agreement with the simulations in Ref. [13 . 
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FIG. 6; Solutions h{t),h(t),C{Z{t),t) for 7 = 0.01, p = 100, 
m = 0.003, and vp = 0.5. 



FIG. 7: Solutions h(t),h(t),C{Z{t),t) for 7 = 0.01, p = 100, 
m = 0.003, and vp = 0.3. 



DISCUSSION 

A crucial point of our analysis is the observation that 
the periodic motion of a planar solidification can only be 
explained, if we go one step beyond the quasi-stationary 



approximation in the expansion ( 25 ) . For a quantitative 



evaluation higher-order terms can be neglected, because 
the definition of the effective expansion parameter m? 
in Eqs. (10) implies m? 5 x 10~^, if we adopt from 



Ref. ^8 the value S — 2 x lO^K/cm of the temperature 
gradient, and from Ref. [T^] the material parameters 
Tm = 1728 K,L = 2350 J/cm3,cr = 3.7 x IQ-^ J/cm^ for 



Nickel, and the interface thickness 2^ = 1.68 x 10 



The fact that the restoring force in Eq. (331 is given 
by iv?h{t) may raise the suspicion that the presence of 
a temperature gradient is an essential ingredient of our 
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FIG. 8: Orbit of the cycle for 7 = 0.01, p = 100, m = 0.003, 
and vp — 0.3. 



0.001 




0.001 



0.01 



0.1 7 



FIG. 9: Neutral lines, enclosing the regions of the Cahn (solid 
line), and of the Mullins-Sekerka (dashed line) instability. The 
vertical line is the projection of the limit cycles in Fig. 9. 



theory. This is only true, however, for a planar geometry 
of the solidification front. In case of a growing spherical 
nucleus the parameter turns out to be proportional 
to the ratio £,/Rc where Rc is the critical radius of the 
droplet. 

Generally, the parameter m determines the period Qq 
of a limit cycle. Close to the threshold at vc the friction 



term in Eq. ( 33 ) can be neglected, so that, due to scaling, 
rig ni. Deep inside the limit-cycle regime accelerations 
are negligible in most parts of the trajectories h{t) in 
Figs. 7 and 10, suggesting to neglect the inertial term 
in Eq. (33 1. Its scaling behavior then implies f2o oc m^, 
in accordance with the statement in Ref. [8] that the 
band width in a pronounced banded structure is inversely 
proportional to the temperature gradient. 

As a final point we note that most phcnomenological 
approaches are based on the assumption of an iV-shaped 
force-velocity relation. This suggests the formation of a 
hysteresis loop, which is considered to represent the limit 
cycle, describing the defect oscillations. In our model the 
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FIG. 10: Solutions h{t), h(t), C{Z{t), t) for 7 = 0.02, p = 100, 
m = 0.003, and vp = 0.5. 



driving force is a convex function of velocity, excluding 
the existence of a hysteresis loop. Instead the necessary 
turnaround of a trajectory at low velocities is provided by 
the inertial term in the oscillator equation, which proves 
the importance of including this term in the equation for 
the interface position. 
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FIG. 11: Neutral lines of the Cahn (solid line), and of the 
Mullins-Sekerka (dashed line) instability. The vertical lines 
are the projections of the limit cycles in Figs. 6 and 10. 



APPENDIX A 

Within our capillary-wave approach the most general 
effective Hamiltonian for the directional solidification of 
a dilute binary alloy reads 



n = / d^x-{dZy 



C-U{z-Z) 



V C -U{z- Z) C-U{z- Z) 
2 



C-U{z-Z) 



(45) 



Here, we have introduced a field C{r,t), which is related 
to the energy density E{r,t) by the equation 



C(.0.|1-£^)S^. (40) 



K L 



Tm 



From this and the equilibrium condition SH/6C = we 
conclude that U{z — Z) obeys the relation 



U{+oo) - U{-oo) 



,47) 

K L J Tm 



which, remembering the relations ([s]) and (|4]), suggests 
to refine our model by assuming 



For the derivation of the model ([l]) - ([T]) the physical 
meanings of the coupling constants i>,k are irrelevant, 
because they will be absorbed into renormalizations of 



the parameters kq and Dq. The only generally important 
constraint on the coupling constants is 



kqk — V > 



(49) 



which ensures stability of the Hamiltonian ( 45 1 



The equations of motion of the generalized model read 

SH 



dtZ 
dtC 



A 



SZ ' 



(50) 



dtC = DV' 



,ji_ sn 

1 5H 



K SC 

where D is the heat diffusion constant. The relation 



r(r,t) = rs- 



sn 

5C 



(51) 



defines a temperature field via a shifted local Legendre 
transform of C(r, i), obeying the condition T{r,t) — Ts 
in thermal equilibrium SH/SC = 0. 

In the limiting case of an infini te h eat conductivity, 
D — oo, and the last of the Eqs. (50 1 is solved, for the 
boundary conditions T{Zp) = Tp, T'{Zp) = S, by the 
static temperature field 



T{z)=Tp + S{z- Zp) 



(52) 



Insertion of the Hamiltonian ( 45 1 and the result ( 52 1 into 



Eq. (51) leads to the relation 



K Civ, t) -U{z- Z) +v C(r, t) - U{z - Z) 



Ts-Tp + S{z - 



(53) 



If the expression for C —U, extracted from Eq. (53 1, 



is inserted into the second of the Eqs. ([50j), one recovers 
the corresponding equation in Eqs. ([5| with the reduced 
diffusion constant 



= I 1 ^]Do. 

KqK^ 



(54) 



In the calculation of the force — SH/SZ, entering the first 



of the Eqs. (50), the relations (48) and (53) can be used 



to eliminate the quantities U, C, which leads to the result 
SH 



SZ 



ad^Z 



(55) 



-\-oo 

-K I dz U'{z - Z) C-U{z- Z) 

-OO 

L '■+°° 
Tm^C 



dz U'{z - Z) Ts-Tp + S{z - Zp) 
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with the renormalized coupling constant 



K = Kq — 



(56) 



The first two terms on the right-hand side of Eq. (55 1 



are identical to the force — SH/ SZ in the first equation 
in Eqs. ([5|. Assuming that J7'(C) is an even function, as 
in case of the model ( 39 ) , the last integral in Eq. ( 55 ) 
reduces to the driving force (|6]). We mention that the 
coupling term cx in the Hamiltonian ([T]) only gives rise 
to a shift in Eqs. (46 1, (47), (48), and to parameter 



renormalizations in Eqs. (54), (56), which all disappear 
in the commonly considered case v = 0. 



APPENDIX B 



Following Ref. [TB], we look for solutions of Eq. (61) 
in form of the expansion 



f{t) = a{t) cos tjj{t) + e''uy{a{t),i!{t)) , (63) 

1^=1,2 

complemented by the constraints 



dtp Ui,{a, ip) simp = 



2tt 



dip u^{a, ip) cos "0 = 0, 



(64) 



which ensure that the u^{a,ip) for u = 1,2 only contain 
higher harmonics in -0. The equations 



In order to derive the differential equation ( 42 ) for the 



amplitude a{t), we start from the oscillator equation (33 ), 
rewritten in the form 



M{h) 



1 1 

M(0) " M{h) 



h . (57) 



where fl, R{h) and M{h) depend parametrically on vp. 
Close to the stability threshold of Eq. (38) is sufficient 



to evaluate all terms in Eq. (57) to leading order of an 
expansion in 



Vp -vc , F'p{vc) = 0. 



(58) 



Since, according to Eq. (38), the linear part of R{h) in h 



is of order e, the leading expression for Eq. ([57j) has the 
general structure 

h + n\vc)h^eXi{vc)h + Z{h,h;vc) ■ (59) 
If we, finally, apply the scaling transformation 

h{t) = ef{t), (60) 
we obtain, up to second order in e, the representation 



(61) 



which is the starting point of the work by Bogoliubov 
and Mitropolsky [16j where in the present case 



?i(/,/) = X,f + X2f'+Y,ff 



(62) 



da 
~di 

dip 
~dt 



i/=l,2 



(65) 



also assumed in Ref. [Tl|, reflect the conditions that the 
amplitude a{t) and the difference 'p{t) — Vlt are slowly 
varying variables. 

We are mainly interested in the functions A^{a) and 



Bv{a), which can be obtained by projecting Eq. (61) 



onto the first harmonics smip and cost/). Then, due to 



Eqs. (64), the contributions u^{a,'p) in Eq. (63) cancel. 



which allows us to look from the beginning for solutions 
of the simplified form 



f{t) = a{t) cos iP{t) . (66) 
Its first and second derivatives are given by 



j{t) = -VLa{t)siniP{t) 

+ E^ [A^cosip — aB^siTLip) , 



(67) 



f{t) = -n"^ a{t) cos iP{t) (68) 
- 2rj E e''(A^sin?A + aBi, cos ip) 



v=1.2 



. dAi „o , 

Ai aBf cos Ip 

da 



dB 

2AiBi + aAi—-^ ] sinip 
da 



From this and the relations (62 1 we, finally, obtain 
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-Bi(a) = - 



1 



gi(/,/))cosV = 



(69) 



AUa) 



in- 1 

Qi(/,/))sinV' = 7^10; 



^2 (a) = - 



27r 
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After scaling back to the variable h{t) via Eq. (60), one 



recovers the result ( 42 ) , where 



a{t) = ea{t) 



(70) 



The coefficients ri, p^, appearing in this equation, depend 
on the parameters 7,p, m, and on the critical velocity vc. 



and can be calculated from Eqs. ( 59 1 , ( 34 1 , and ( 62 1 . For 
the choice 7 = 0.01, p = 100, m = 0.003, and vc = 0.5214 
one finds n = 1.5575 • 10"^ and pa = 9.885 • 10"^ 
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